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ABSTRACT 

We make use of Spitzer imaging between 4 and 16 /im and near-infrared data at 
2.2 /zm to investigate the nature and distribution of the mid-infrared emission in a 
sample of early-type galaxies in the Virgo cluster. These data allow us to conclude, 
with some confidence, that the emission at 16 /im in passive ETGs is stellar in ori- 
gin, consistent with previous work concluding that the excess mid-infrared emission 
comes from the dusty envelopes around evolved AGB stars. There is little evidence 
for the mid- infrared emission of an unresolved central component, as might arise in 
the presence of a dusty torus associated with a low-luminosity AGN. 

We nonetheless find that the 16 /im emission is more centrally peaked than the 
near- infrared emission, implying a radial stellar population gradient. By comparing 
with independent evidence from studies at optical wavelengths, we conclude that a 
metallicity that falls with increasing radius is the principal driver of the observed 
gradient. 

We also plot the mid-infrared colour-magnitude diagram and combine with similar 
work on the Coma cluster to define the colour-magnitude relation for absolute K-band 
magnitudes from -26 to -19. Because a correlation between mass and age would produce 
a relation with a gradient in the opposite sense to that observed, we conclude that the 
relation reflects the fact that passive ETGs of lower mass also have a lower average 
metallicity. The colour-magnitude relation is thus driven by metallicity effects. 

In contrast to what is found in Coma, we do not find any objects with anomalously 
bright 16 /urn emission relative to the colour-magnitude relation. Although there is 
little overlap in the mass ranges probed in the two clusters, this may suggest that 
observable "rejuvenation" episodes are limited to intermediate mass objects. 

Key words: galaxies: elliptical and lenticular cD, galaxies:evolution, galax- 
ies:clusters:general, galaxies:photometry, infrared:galaxies. 



1 INTRODUCTION 

The mid-infrared emission of even the most passive early- 
type galaxies (ETGs) shows an excess of emission longward 
of ~ 8 (im over that which is expected from purely photo- 
spheric emission. This excess emission was detected by Im- 
pey et al. (1986) using ground-based observations and was 
subsequently studied using ISO (e.g. Bressan et al. 2001). 
The excess has now been observed in the central regions of 
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bright ETGs in the Virgo cluster (Bressan et al. 2006) us- 
ing Spitzer-IRS. Although some ETGs show a mid-infrared 
excess that is evidently caused by warm dust in a star form- 
ing inter-stellar medium (Panuzzo et al. 2007) the excess in 
passive objects can be explained as the integrated emission 
from the hot dust in the envelopes of evolved AGB stars. 
However, hot dust emission could also arise from the dusty 
torus around a central active galactic nucleus (AGN) of low 
luminosity. 

Although there is evidence that the emission at 16 /im 
is extended in Coma cluster galaxies (Clemens et al. 2009a) 
the weight of evidence for this is far from conclusive. In the 
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present article we discuss Spitzer-IRS peakup imaging at 
16 (im of ETGs in the Virgo cluster that permit the investi- 
gation of the spatial distribution of the excess mid-infrared 
emission. In order to isolate the effects of the mid-infrared 
excess caused by AGB stars we also use IRAC images at 
4.5 and 8.0 /tm because these are relatively unaffected by 
the excess and sample the purely photospheric component 
at the longest possible wavelength. We will argue that the 
excess is stellar in origin, as has been previously suggested 
(Bressan et al. 1998, Athey et al. 2002, Xilouris et al. 2003). 

The stellar origin of the excess mid-infrared emission 
makes it a particularly useful age/metallicity diagnostic be- 
cause the contribution of the mid-infrared excess to the inte- 
grated spectrum varies with age and metallicity in a different 
way to optical diagnostics (Bressan, Granato & Silva 1998). 
As a stellar population gets younger and/or the metallic- 
ity increases the mid-infrared excess increases. However, the 
optical H/3 index becomes larger (optical colours bluer) as 
the age decreases and/or the metallicity decreases. The op- 
posite behavior of optical indices (colours, narrow band in- 
dices and spectral shape) and mid-IR excess with respect to 
age and metallicity variations, means that the mid-infrared 
spectral region, in combination with optical data, can dis- 
entangle age-metallicity effects in early-type galaxies. Mid- 
infrared observations, then, can provide without the need 
of accurate modeling the explanation of the origin of the 
colour-magnitude relation in clusters of galaxies. This can 
be understood by considering the two following extremes. 

• If galaxies are coeval and the bluer colours toward the 
less luminous sources are due to lower average metallicity 
(e.g. Arimoto & Yoshii, 1986) we expect the mid-infrared 
excess to become relatively less pronounced at decreasing 
luminosity. 

• If the galaxies have the same average metallicity and 
the bluer colours toward less luminous sources are due to a 
younger age (the alternative explanation allowed in principle 
by optical diagnostics) we expect the mid-infrared excess to 
become relatively more pronounced at decreasing luminos- 
ity. 

Although most studies tend to find that both age and 
metallicity vary as a function of galaxy luminosity or veloc- 
ity dispersion (e.g. Nelan et al. 2005, Thomas et al. 2005, 
Clemens et al. 2006, 2009b) it is not clear which is the main 
driver of the colour-magnitude relation in the cluster envi- 
ronment. Here we attempt to address this question using a 
multi-wavelength appraoch. 

In Clemens et al. (2009a) we used Spitzer-IRS peakup 
imaging at 16 /im to populate the luminous part of the mid- 
infrared colour magnitude diagram for the Coma cluster. In 
Coma, galaxies with K-band magnitudes only as faint as -22 
could be studied. In the present paper we repeat the analysis 
for fainter objects in the Virgo cluster in order to extend the 
colour-magnitude relation to fainter magnitudes. 



2 DATA ANALYSIS 

The aperture photometry for both the IRAC and IRS- 
peakup images is based on post basic calibration data 
(PBCD). 



The 16 /xm, blue peakup images are background domi- 
nated and the background often shows fluctuations on spa- 
tial scales similar to the source sizes. Background subtrac- 
tion techniques that attempt to model the background as 
the smoothly varying component of an image thus either in- 
clude extended emission from the galaxy as background or 
do not correct for the background fluctuations. This tends 
to cause an over-subtraction of the background in the vicin- 
ity of the galaxy. We therefore define the background level, 
at all wavelengths, to be the median value in an elliptical 
annulus, of width 1.5 pixels, centred on the galaxy, with 
semi-major axis and position angle as given by the mean 
optical diameter catalogued in "Hyperleda'Q. Thus we sam- 
ple the background at twice the radius of the galaxy. This is 
actually the background estimation technique recommended 
for IRAC images of extended sourcefl 

Confusing sources, normally in the background, can be 
problematic, especially at 16 /jm. For meaningful aperture 
aperture photometry they need to be removed. In some cases 
the contaminating sources were too close or bright to at- 
tempt removal (e.g. vcc 0951). In some cases the 16 /im 
detections themselves were not at the catalogued position 
and so were taken as "suspect" background sources and dis- 
carded (e.g. NGC 4366). For the remaining objects that 
showed nearby contaminating sources we applied an algo- 
rithm to remove them which analyzed the pixel distribution 
in a series of annuli with axis ratio and position angle equal 
to that of the galaxies. Pixels that deviated by more than a 
fixed number of standard deviations from the median value 
were replaced by the median value. The threshold for source 
removal was varied between two and four standard devia- 
tions and was sometimes applied iteratively to achieve the 
best results. In all cases the results were carefully checked 
by eye. 

Integrated fluxes were taken as the total flux within the 
background annulus. We stress that these fluxes are mea- 
sured only for the purposes of determining the mean colour 
of the galaxies. Given the way in which the background is 
determined, they are not a good measure of the true inte- 
grated fluxes of the objects. The colours of the objects, how- 
ever, should be robustly determined. When we come to plot 
the colour-magnitude diagram we will use the total K-band 
magnitude of the galaxies as given in the 2MASS catalogue, 
rather than our measurement from the 2MASS K-band im- 
age. 

Radial colour profiles were determined using a series 
of elliptical annuli with axis ratios and position angles as 
given in "Hyperleda" . Before fluxes were measured the two 
images were first convolved to identical resolutions. This 
was achieved by convolving each image with the measured 
PSF of the other. For example, when considering the K-[16] 
colour the 2MASS K-band image was convolved with the 
Spitzer blue peakup PSF and the 16 /im image was con- 
volved with the 2MASS K-band PSF. Instrumental resolu- 
tions (FWHM) are 2'.'6, l'.'7, l'.'9 and 3'.'6 for 2MASS, IRAC 
4.5 and 8 and 16 respectively. 

The error on the derived fluxes was taken as the quadra- 
ture sum of the rms variation within the aperture or annulus 
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and the calibration error, which was taken as 5% for the blue 
peakup data and 10% for the IRAC data (see the respective 
instrument handbooks). 

All aperture photometry was performed using a custom 
pipeline written in idl and use was made of the IDL As- 
tronomy User's Library (Landsman 1995). 

3 RESULTS 

3.1 Contribution from a central point source in 
the mid-infrared? 

Most ETGs in our sample show negative radial gradients in 
the K-[16] colour, indicating that the mid-infrared 16 /an 
emission is more centrally concentrated than the K-band 
light. The IRS-spectroscopy of Bressan et al. (2006) has 
already shown this emission to be extended. However, at 
16 )im there are roughly similar contributions from the pho- 
tospheric and dusty envelope components (See figure 1 of 
Clemens et al. 2010) so it remains possible that only the 
stellar component is extended, while the dusty, excess, com- 
ponent is unresolved. This might be expected if these galax- 
ies contained a dusty torus at their centre, associated with 
a low luminosity AGN. 

Using HST images, Carollo et al. (1997) interpreted the 
central dip in the optical surface brightness of NGC 4406 
within 0'.'4 as evidence of a nuclear dust ring. NGC 4636 is 
classified as a LINER/Sy 3 and has radio jets at 1.4 and 
4.8 GHz (Stanger & Warwick, 1986). Both of these objects, 
then, show evidence of AGN activity. So is there any emis- 
sion from a dusty torus at 16 /im? 

Although Clemens et al. (2009a) found no evidence for 
central point sources at 16 /tni for the galaxies of Coma, 
based on effective radii estimates, the spatial resolution was 
insufficient to investigate this possibility in detail. 

In order to test this possibility we created model 16 /im 
images based on the K-band images for the four large galax- 
ies in our sample. We first convolved the K-band images 
to the same resolution as the observed 16 /im images, and 
then added a central point source in the form of the 16 /jm 
psf. We then measured the radial surface brightness profiles 
for point sources of various strengths and compared with 
the observed 16 /im profiles. The results are illustrated in 

Kg. m 

We normalized the models in two different ways. Firstly, 
to the observed surface brightness at a radius of 15", so as to 
be far from any significant PSF effects, and secondly so that 
the observed central surface brightness was equal to that of 
the model with a 60% contribution from a point source. The 
choice of 60% is based on model fits to Spitzer-IRS profiles 
of passive early-type galaxies that show that approximately 
this fraction of the 16 /im emission comes from the dusty 
circumstellar envelopes of AGB stars (Bressan et al. 2006) . If 
the colour gradient is due to a population gradient involving 
AGB stars then this would be a sensible upper limit to the 
extent to which the colour of the nucleus might differ from 
that at larger radii. 

Fig. [2] shows that the observed profiles tend to have 
shallower gradients than one would expect if a central point 
source were the cause of the observed colour gradients. 
Nonetheless, the limited signal-to-noise ratio and spatial res- 
olution of the 16 /xm images means that a central point 



source is permitted at some level. The most interesting case 
is NGC 4636 because this contains a radio loud AGN. The 
observed profile in Fig. [5] is clearly much shallower than 
would be expected if there were a significant contribution 
from a central point source. At least for this galaxy, there is 
no evidence of a significant contribution to the mid-infrared 
emission from warm dust in a central torus. The observed 
colour gradient is therefore due to a stellar population gradi- 
ent. If this is so for an ETG with a known AGN, there seems 
little evidence that the colour gradients in the other galaxies 
are due to anything else. A similar conclusion was actually 
reached for the mid-infrared emission of M 87, where the 
only excess over purely stellar emission is synchrotron emis- 
sion from the radio loud nucleus (Buson et al. 2009, Baes et 
al. 2010). 

3.2 Radial colour gradients 

Fig. [3]shows the radial mid-infrared colour profiles, K-[4.5], 
K-[8] and K-[16]. The signal-to-noise ratio of the 16 /im im- 
ages limits the number of galaxies for which the K-[16] colour 
profile can be plotted, and in general, profiles are truncated 
before noise becomes dominant. In addition, any objects 
that had nearby confusing sources have been excluded from 
the analysis. The 4 brightest galaxies are separated in the 
figure and have profiles that are well measured within a ra- 
dius of ~ 1'. We place more emphasis on these objects in 
the following discussion as the profiles are the most reliable. 

For the 4 brightest objects, and for the majority of the 
fainter galaxies, at least at radii less than 10", the K-[16] 
colour shows a negative radial gradient that is steeper than 
that of either the K-[4.5] or K-[8] gradients. We also see 
that (with the exception of NGC 4516) the K-[4.5] and K-[8] 
profiles are very nearly parallel. In several cases, including 
the brightest galaxies, these profiles are almost flat. 

With reference to the model fits to mid-infrared Spitzer- 
IRS spectra (Bressan et al. 2006) we expect both the 4.5 and 
8 /mi micron fluxes to be dominated by emission from stel- 
lar photospheres. However, at 16 /im this changes, with ap- 
proximately 60% of the emission coming from the hot dusty 
envelopes around evolved AGB stars. This offers a very natu- 
ral explanation for the observed difference between the mid- 
infrared colour profiles, suggesting that the radial distribu- 
tion of dusty AGB stars is more centrally concentrated than 
that of the general stellar population. 

3.3 Colour— magnitude diagram 

Table[T]gives the mid-infrared colours of our sample galaxies 
determined within the optical radius, as described in sec- 
tion 

In Fig. [4] we show the mid-infrared colour-magnitude 
diagram for those Virgo galaxies for which the K-[16] colour 
could be reliably determined. In light grey we reproduce 
the values obtained for the Coma cluster by Clemens et al. 
(2009a). The galaxies observed in the Virgo cluster allow 
this diagram to be extended to fainter magnitudes than was 
possible for the Coma cluster alone, with objects down to 
an absolute K-band magnitude of -19.23. 

The colour magnitude relation of the Coma cluster 
showed a trend for bluer K-[16] colours towards fainter mag- 
nitudes. The addition of the fainter Virgo galaxies shows 
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Figure 1. Example 16 fim images. The annuli used for the determination of the background are shown. Background sources that can 
be seen in some images were removed before the photometric measurements. 
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Figure 4. MIR colour-magnitude diagram. The colours are de- 
fined as those within 80% of the optical radius given in "Hyper- 
leda" . K-band fluxes were taken from the 2MASS All Sky Survey. 



that this trend continues to fainter magnitudes. There ap- 
pears to be a slight offset between the Coma and Virgo data. 
Although this could be caused by an error in the relative 
distance of the two clusters (taken to be 97 and 17 Mpc re- 
spectively) it may also be due to the fact that the colours 
were determined within larger apertures for the Coma clus- 
ter. The sense of the offset (Virgo objects appearing with 
slightly higher than expected K-[16]) is consistent with the 
colour gradients described in section [3721 

We also note that the 4 bright Virgo galaxies lie on 
the low boundary of values of K-[16] found in the Coma 
cluster. There is no obvious data analysis effect that could 
have caused such a bias. 



4 DISCUSSION 

Both age and metallicity can effect the strength of the in- 
tegrated 16 fim emission from the dusty envelopes of the 
AGB population. Younger stellar populations, for example, 
should have a larger contribution from dusty AGB stellar en- 
velopes. However, in order to separate the effects of age and 
metallicity we will compare our results with those of studies 
at optical wavelengths. As discussed in section [TJ the differ- 
ent responses of observables in the optical and infrared, to 
changes in metallicity and age, can be used to resolve the 
degeneracy between these two parameters. 



4.1 Radial colour gradients 



As discussed in section 13.11 it seems likely that the radial 
gradients in the K-[16] colour are due to a population gradi- 
ent of AGB stars relative to the general stellar population. 
The redder K-[16] colours toward the centre of the galaxies 
could be due to younger ages, higher metallicities or both. 

Clemens et al. (2009b) have used optical line strength 
indices of a sample of 14,000 ETGs from the Sloan Dig- 
ital Sky Survey (SDSS) to determine age and metallicity 
as a function of mass, environment and also galactic ra- 
dius. They find negative radial metallicity gradients for all 
masses. Although they also find positive radial age gradients 
for massive galaxies this trend is less apparent at interme- 
diate masses and may even invert for low mass galaxies. 
The 4 bright galaxies in our Virgo cluster sample have cen- 
tral velocity dispersions in the range 170-250 kms" 1 . In this 
case little or no age gradient is expected, at least on a sta- 
tistical basis. Although it must be borne in mind that the 
results based on the SDSS are statistical in nature, and not 
necessarily applicable to any individual object, these results 
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Figure 2. Comparison of model 16 fira surface brightness profiles (solid black lines) with the observed profiles for the four largest galaxies 
in the sample (dashed lines). For the red line the models and observed profile have been normalized at a radius of 15", well beyond any 
significant psf effects. The blue line, instead, is normalized to a central point source contribution of 60%, approximately the contribution 
of dusty AGB star envelopes as measured in Spitzer-IRS spectra. The model profiles assume a 16 fim brightness similar to that in the 
K-band with an additional point source at the centre. The contribution of the point source rises from zero (bottom line) in increments 
of 10%. Surface brightness units are arbitrary. 



would favour metallicity effects over age effects as the cause 
of the mid-infrared colour gradients. 

Negative metallicity gradients and negligible age gradi- 
ents have, in fact, been found by other studies of ETGs both 
in the field (Annibali et al, 2007) and cluster (Rawle et al. 
2008, Mehlert et al. 2003, Peletier et al. 1999) using optical 
line strength indices. 



4.2 Colour-magnitude relation 

Since the work of Bower, Lucey & Ellis (1992) the colour- 
magnitude relation (its slope and narrowness) has been open 
to two interpretations. One is that ETGs are all very old and 
the slope of the colour-magnitude relation reflects a trend 
for lower mass galaxies to be less metal-rich. The narrowness 
results from a small dispersion in the formation epoch with 
respect to the age of the system. In the second interpreta- 
tion, if the luminosity weighted "average age" is allowed to 
vary, the only way to maintain a small dispersion is to in- 
voke a synchronization mechanism, in such a way that lower 
mass ETGs appear younger. 



For many years, the metallicity trend had been favoured 
partly because of its simplicity and partly because of lack of 
evidence for a relation between mass and age (e.g Kodama 
et al. 1998). However, recently, numerous studies have found 
evidence for a so-called "downsizing" effect (e.g. Cimatti, 
Daddi & Rezini, 2006), which is nothing less than such a 
relation, in the sense that star formation is more prolonged 
in galaxies of lower mass. The evidence of the presence of a 
mechanism that links the mean age of an object to its mass 
restores some of the original doubt concerning the cause of 
the narrow colour-magnitude relation. 

Studies at optical wavelengths, based on line strength 
indices, find that less massive galaxies tend to be younger 
and less metal rich (Bressan et al. 1996, Kuntschner et al. 
2001, Thomas et al. 2005, Nelan et al. 2005, Clemens et al. 
2006, Smith et al. 2007), both in the field and in the cluster 
environment. If the age trend holds also for the galaxies in 
our sample (both in Virgo and in Coma) then for a fixed 
metallicity, we would expect lower mass galaxies to have 
more 16 /zm emission per K-band flux relative to the more 
massive objects, because of the greater importance of the 
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Figure 3. Mid-infrared radial colour profiles. Red: K s — [16], grcen:K a — [8], blue: K s — [4.5]. 



dusty AGB phase. That is, low luminosity objects should 
have higher values of K-[16] if age were the dominant effect. 
The observed colour-magnitude relation, however, shows the 
opposite trend, with lower mass systems having lower K- 
[16] colours. Metallicity, then, drives the mid-infrared colour- 
magnitude relation such that stellar populations with lower 
metallicity emit less 16 /jm flux per unit stellar mass than 
their more metal-rich and more massive counterparts. In 
addition, population synthesis modelling of the K-[16] vs 



V-K colour plane shows that the optical colour-magnitude 
relation is also driven by a dominant metallicity effect (see 
Fig. 6 of Clemens et al. 2009a). Thus, Spitzer shows that the 
cluster colour-magnitude relation is a sequence of metallicity 
rather than age. 

Therefore, although a relation between mass and age 
could exist for our sample galaxies this effect is dominated by 
a mass-metallicity trend that defines the colour-magnitude 
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Table 1. Mean mid-infrared clours. Missing values imply either a non-detection in one band or that the source was confused. The errors 
in this table include the systematic calibration errors. These are not included in the error bars of Fig. [4] 



relation. Thus, "downsizing", if present, is of secondary im- 
portance in passive ETGs. 

The colour-magnitude diagram of the Coma cluster 
shows several objects that have high 16 /J,m fluxes relative 
to the colour-magnitude relation. None of the ETGs in the 
present Virgo sample show this. Although the statistics are 
rather poorer for the Virgo sample and in Virgo we do not 
sample the luminosity range where these anomalous colours 
are seen in Coma, this may nonetheless indicate that the re- 
juvenation events that can be detected in the mid-infrared 
due to the increased 16 /tm emission, are limited to inter- 
mediate mass objects. 



5 CONCLUSIONS 

We have observed a sample of low-luminosity ETGs in the 
Virgo cluster using the blue peakup detector of Spitzer-IRS. 
Our aperture photometry of the spatially resolved images 
can be used to draw the following conclusions. 

(i) The 16 /xm emission that is in excess of that from 
stellar photospheres is extended. This emission is therefore 



consistent with an association with the dusty circumstellar 
envelopes around evolved AGB stars, rather than from a 
dusty torus around a low luminosity AGN. This is a more 
robust conclusion than earlier evidence based on Spitzer-IRS 
spectroscopy (Bressan et al. 2006). 

(ii) Mid-infrared colour gradients tend to show that the 
16 /Urn emission is more centrally concentrated than the K- 
band light. Much shallower colour gradients involving bands 
dominated by photospheric emission, K-[4.5] and K-[8], im- 
ply a population gradient in the AGB stellar population. 
Comparison with independent studies at optical wavelengths 
implies that this is a metallicity effect. 

(iii) The mid-infrared colour magnitude relation, that we 
define over 8 magnitudes by combining with similar data 
for the Coma cluster, is driven by metallicity. Less massive 
objects are less metal rich. 

(iv) If there is a relation between galaxy mass and age in 
these clusters the effect is masked by the dominant metallic- 
ity effect. Thus, the mid-infrared colour-magnitude relation, 
together with optical studies, shows no evidence for so-called 
"downsizing" in the cluster environment. 

(v) The lack of objects with K-[16] colours that place 
them well-above the colour-magnitude relation of passive 
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ETGs in the present Virgo sample, and for the most mas- 
sive objects in the Coma sample of Clemens et al. (2009a), 
may be evidence that observable rejuvenation episodes are 
limited to intermediate mass objects. 
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